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3 Abstract 

I— 9 

New vectorlike fermions that mix with the third generation can significantly affect the r and b 

m 

Yukawa couplings. Consistent with precision electroweak measurements, the width of the Higgs 
boson to rr, bb can be reduced by Oil) with respect to the Standard Model values. In the case 

i 

Q_i of the b quark, a reduced width would result in an enhanced branching ratio for other final states, 

<D 

such as 77. New leptons can also substantially modify the Higgs boson branching ratio to photons 
through radiative effects, while new quarks can contribute to gg fusion. The combined effect can 

> 

be as much as a factor of two on the branching ratio to 77. The new quarks and leptons could be 

light, which would allow discovery at the LHC. In the case of significant suppression of h —> tt, 

^1 searches for new leptons decaying to r-rich final states, perhaps in association with Higgs bosons 

O 

(S) are motivated. 

•i-H 
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I. INTRODUCTION 



Following the observation of a "Higgs-like" state near 125 GeV [HE], it remains to precisely 
determine whether the couplings are in fact those expected from a Standard Model (SM) 
Higgs boson. Branching ratios to a variety of final states are non-trivial at this mass 
value - many measurements can be made and compared to their corresponding SM pre- 
dictions. Deviations from these measurements could potentially indicate the presence of new 
physics. 

If the Standard Model is a good effective theory near the weak scale, a modification to 
the width to fermions T(h — > ff) is realized through the presence of the dimension-six 
operator 

O h s = (f D hf c )(rfh), (1) 

where fu represents an SU(2) doublet, and f c is the right-handed partner. When combined 
with a Standard Model-like Yukawa coupling yjfohf , the mass and effective Higgs Yukawa 
coupling of the / are given by 

m/ = yjv + c h sv 3 , (2) 
yf = y° f + 3c h3 v 2 , (3) 

where v = 174 GeV and represents the coefficient of O^- The mismatch between Eqs. (|2]) 
and (|3j) indicates the possibility for a discrepancy between the observed fermion mass and 
the Yukawa coupling. Of particular interest are modifications for f = r,b. These couplings 
are small enough that it is plausible for them to be affected by integrating out new physics 
near the weak scale, but not so small that it is hopeless to measure them in the near future. 
In this paper, we concentrate on the modification of these two couplings. 

An interesting secondary effect arises for f = b. Since for = 125.5 GeV the b makes up 
a large fraction of the total width, a suppression (or enhancement) of T(h — > bb) will affect 
the branching ratio (BR) to all other states. Defining Rh^bb = r(/i — > bb)/Y SM {h — > bb), in 
Fig. [l] we illustrate the effect of a modification of the b width on the BR to bb (solid) and 
to all other states (dashed). For example, a 40% reduction in T(h — > bb) can give rise to a 
~ 30% enhancement in the branching ratios of all other channels. 

Previously, realizations of this operator have been considered in the context of two Higgs 
doublet models (e.g., |4j). Here, we realize the operator of Eq. ([!]) by integrating out 
new vectorlike fermions. We examine how large its effect can be consistent with precision 
electroweak constraints and discuss its implications for Higgs phenomenology. We highlight 
how these new vectorlike fermions could be of possible interest for two trends realized in 
the current data: the branching ratio to 77 seems somewhat higher than expected [H [2], 
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FIG. 1. The effect of a modification of the width of a 125.5 GeV Higgs boson to b quarks (plotted 
as a function of Rh-^bb = — > bb) /T sAi(h — > bb)) on the branching ratio to b quarks (solid) and 
branching ratios to all other states (dashed). The Standard Model widths are taken from [3]. 



while the branching ratio the rr is lower than expected [5], see, e.g., [6]. In contrast to 
two Higgs doublet models where, for example, the observed Higgs has minimal overlap with 
the down- type Higgs boson [7j, here the r Yukawa may be modified independently of the 
b Yukawa. Vectorlike leptons allow the generation of the operator of Eq. Q with / = r 
without the corresponding operator for f = b. This may be of interest given the hint of a 
signal in bb at the Tevatron [8]. 

The remainder of the paper is organized as follows. In section [TTJ, we describe a model 
involving vectorlike fermions that can give rise to the operator of Eq. Q. We focus on two 
cases in which this leads to the modification of the Higgs boson coupling to r and b. In 



section III we discuss other effective field theory operators generated in this model that can 



be used to place constraints on the new physics. These constraints will determine the size of 



the Higgs coupling modifications it is possible to achieve. In section [IV] we discuss regions 
of parameter space in which these new states may also affect the effective Higgs coupling to 
massless gauge bosons, notably highlighting how the vectorlike leptons can sizably enhance 
h — > 77. Results related to this point have recently appeared in [9HT2]. We briefly explore 
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discovery possibilities for the new states at the LHC in section |Vj and comment on possible 
UV completions in section VI Finally, our conclusions are presented in section VII 



II. THE MODEL 



One way to realize the operator of Eq. (jl]) is by mixing the Standard Model fermions with 
new heavy fermions. As an example, we focus first on leptons. These new leptons have 
additional sources of mass and, upon mixing, the r inherits some of this mass. 



We write the "Standard Model" tau lepton doublet as 




(4) 



and the corresponding SU (2) L singlet field as e c . We then augment the Standard Model by 
a vectorlike pair of SU (2) L doublets 




and L 



N 



(5) 



and a vectorlike pair of SU(2) E singlets, E c and E c . Both for simplicity and motivated by 
flavor constraints, we ignore mixing with the first two generations of leptons. The presence of 
these fields allows several new Yukawa couplings. The mass terms and interactions are 

- C 3 y° T ie c h + y E lE c h + y L Le c h + y LE LE c h + y LE LE c h ] + M E E C E C + M L LL + h.c. (6) 

We have rotated away possible terms of the form [iitL and fi e e c E c and labelled Yukawa 
couplings by the exotic fermion(s) present in the interaction. When the Higgs field is set 
to its vacuum expectation value, this leads to mass terms for the charged leptons of the 
form 



-A, 



e c E c L + )M 



E° 

\ L -> 



f m o 



+ h.c, where Ai 



\ 



(7) 



y L v 
y E v M E y LE v 
y y LE v M L J 

with v = 174 GeV. M. is diagonalized to yield three charged Dirac fermions - the r plus two 
exotic, charged leptons denoted 
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with i = 1, 2. For Mb 



d = diag(mi, 
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(9) 



The spectrum also contains the massless Standard Model neutrino and a massive neutral 
Dirac fermion (consisting of N and N) with mass Ml- 

An analogous model can be written down to modify the effective Higgs Yukawa coupling of 
the b by making the replacements 

(-) (-) (-) (-) 
I -> q,e c ->■ b c , L^Q , £^L> . (10) 

In this case, we denote the bottom-like quarks as Bi and B 2 respectively, and the top-like 
quark of mass Mq as T. We neglect mixing in the top sector (e.g. via couplings of the 
form Qt ^) for simplicity^] Since we are focused on the possible effects for Higgs physics, 
we assume some alignment that allows us to couple to the third generation exclusively, 
and do not explore flavor models explicitly. The maximal mixing angles between the third 
generation and the new heavy fermions are relatively small for the benchmarks we will 
consider (sin 2 9 < 5 x 10~ 3 ), and so given this assumption the model should be safe from 
flavor constraints. 



III. EFFECTIVE THEORY CONSIDERATIONS 



Integrating out the heavy lepton fields generates a contribution to the effective operator of 
Eq. ([I]) as desired (see Fig. [2]). The leading result is 

Ch3 = . 11 

h M E M L 1 ; 

which can take either sign, allowing it to suppress or enhance the r Yukawa coupling relative 
to the Standard Model value. We will concentrate on the suppression of Yukawa couplings, 
as is presently slightly preferred by the data. Notably, this contribution to O h 3 is not 
proportional to the SM Yukawa coupling, y° ~ 1CT 2 , so can potentially compete with it in 



1 Unless the scale of new physics is very low, mixing with the SM top is unlikely to substantially modify the 
top Yukawa. In any case, this could be effectively absorbed into a modification of the effective coupling 

of the Higgs boson to gluons and a (likely modest) modification of the Higgs boson coupling to photons. 

(-) „ .... (") 

Introduction of JJ fields and their mixings with the Q could further modify the T parameter. 
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FIG. 2. An effective r Yukawa coupling. 



spite of the v 2 /M 2 suppression of Eq. ([2]). Operators of yet higher dimension that include 
novel couplings, e.g. jjle, can be numerically significant (and are included below), but the 
above equation serves as a useful guide to the expected size of the effect 



A. Constraints from the Z-pole 



The above indicates modifications of the r and b Yukawa couplings are possible, but the 
magnitude of the effect clearly depends on allowed size of the Yukawa couplings with the 
exotics. In the case of the leptonic model, the E c mixing with £ (and L mixing with e c ) 
modifies the couplings of the r lepton to gauge bosons, which can provide constraints on 
these couplings. 

We discuss these modifications in an effective theory language p3J E] where the expressions 
in terms of the mass matrix are exceedingly simple. This will serve as an important guide to 
the region of parameter space where large deviations in T(h — > ff) are possible, consistent 
with known experimental constraints. 

The new physics generates operators of the form 

O ht = i(tfD l Ji)(£ T 'y li l r ), (12) 
O'u = i^D^hXl^Q, (13) 
O he = i{tfD il h){f 1 ^T). (14) 

Here r corresponds to the SU(2) singlet part of the r lepton. An effect of these operators is 

2 We also expect additional contributions to c^3 due to modifications of wave-function renormalization of 
the fermions, but these contributions are proportional to and hence negligible. 
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FIG. 3. Diagrams that give the leading contributions to the effective operators of Eq. (12)-(14), 
and hence modification of the Z couplings to the r. 

to shift the couplings to the Z-boson by [I3l [H] 

8g£ = v 2 {c M - c u ) SgX = v 2 (c' M - c M ) (15) 
5g£ = v 2 (c he - CM - c' M ) 8gX = -v 2 (c he + c M + c' M ). (16) 

Such departures in the gauge couplings are constrained by measurements at the Z-pole, in 
particular R T = r(hadrons)/r(rr) as well as the asymmetries A T and A^'^ . The model 
that modifies the b quark Yukawa coupling generates similar operators with £ T — > g& r — > b. 
Note Rb = r(66)/r(hadrons), the reciprocal of the definition for the analogous R T . 

Our new physics generates Om and O'm with identical coefficients so that coupling of the 
neutrino (top quark) is unchanged. Consequently, the shifts in the r couplings simplify: 

8gf = v 2 (c he -2c M ), 8gX = -v 2 (c he + 2c M ). (17) 



Constraints on these values will help determine the maximum size of the effect of Eq. ([11]). 
The diagrams in Fig. [3] generate Ch e and c«. To leading order, the shifts to the r vector and 
axial couplings are given by 



.4 
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2 




Ml 



« = (19) 

Higher order operators (e.g. operators of the form Ohi(H Jf H) n ) can involve ijle and y^E - if 
these couplings are large, their contribution can be relevant. As we perform exact numerical 
diagonalization of the relevant mass matrices, these effects are included in our discussion 
below. Note that while Sg^ has fixed sign, 5g^ can take on either sign (or be tuned small). 



The constraints on Eqs. (18) and (19) from data limit Vl/M l and Ve/Me, and hence will 
limit the size of c^s. With the replacements, L — )■ Q and E — >■ D, the results of this section 
trivially translate to vectorlike quarks. 
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B. Modification of the Tau Yukawa Coupling 



We now turn towards a quantitative discussion of how measurements of r leptons at the 
Z pole constrain the lepton model. We will then be prepared to discuss the size of the 
modifications to the r Yukawa coupling achievable subject to these constraints. 

To leading order in the couplings we have 

Rr oc (g 2 v + gl)- 1 5R T = i^ M (0.3 Sg? + 4.0 Sgf), 

A - K Ifik 5A - = 0.29 5^-3.9 5^, (20) 

A°> T _ 3 ASM A =k A4°' T - 3 ASMS, a 

Experimental results for these quantities, as well as their SM predictions (using the value of 
sm 2 9w found from fitting the entire suite of PEW data) are given by [13] : 

R ex P = 20.764 ±0.045 R S T M = 20.789 ± 0.011, 

A ex P = 0.1439 ±0.0043 Af M = 0.1475 ± 0.0010, (21) 
(Af^yxp = 0.0188 ±0.0017 Af% >T = 0.01633 ± 0.00021. 

The ability to make 5gy small means that, in much of the parameter space, the strongest 
constraint comes from R T . 

For fixed (M L , M E ,y LE ,y LE ), the constraints on bgX and 5g^ can be visualized in terms 



of elliptical Ax 2 contours in the {vl^Ve) plane. Meanwhile, we can see from Eq. (11) that 
lines of constant Rh^rr = r(/i — > TT)/Tsyi(h — > tt) will be approximate hyperbolae in the 
same plane. An example of these curves is shown in Fig. [4] for M E = M E = 350 GeV and 
Vle — Vle — 1- Note the Ax 2 shown is measured with respect to the global minimum in the 
(5g^,Sg^) plane. This minimum has 5g^ < 0, which cannot be achieved in this model, i.e. 
no point with A% 2 = appears on this plot. The largest deviation in R] l ^ tTT consistent with 
requiring that the values of 8g^ and 5g^ give a particular A% 2 can be determined by finding 
the hyperbola of greatest deviation that intersects the appropriate Ax 2 ellipse. For instance, 
in the case of the reference point chosen in Fig. |4| one can achieve Rh^tt ~ 0.7 consistent 
with a (SgV , Sg^) fit satisfying Ax 2 < 5.99. Note that if one chose to allow Ax 2 measured 
instead relative to the Standard Model less than 5.99, it would not significantly alter these 
results, although one could reach slightly smaller values of Rk-^tt ( e -g- Rh-^tt ~ 0.6 for the 
reference point). 

Above, we have taken y EE ~ 0(1) in order to achieve an appreciable affect on Rh^ TT . Even 
larger deviations in Rh^ TT m ay be achieved by increasing y LE . To a good approximation, 
the maximal effect on c/j3 is proportional to this coupling. Taking y LE = y LE = 2 allows 
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FIG. 4. Contours (solid) of Ax 2 = 2.30, 4.61, and 5.99 relative to the best-fit value in the (6g^, Sg^) 
plane as a function of i/l and i/e for M = Me = Me = 350 GeV and i/le = Vle = 1- These values 
of Ax 2 correspond to 68.27%, 90% and 95% regions for 2 parameters. For reference, the Standard 
Model (SgY = 0, 5g^ = or equivalently yi = 0, ue = 0) gives Ax| M = 2.09. Also shown are 
contours (dashed) of Rh^rr = T(/i — > tt)/Tsm(^ — > tt). As discussed in the text, larger values of 
Vle would allow for larger deviations in R^ TT . 



R h ^ TT = 0.45. Couplings this large can produce tension with the isospin breaking parameter 
T that can be ameliorated by going to large M (of order a TeV). We elaborate on this issue 
in section HI D We remain agnostic as to the new physics that would be required at low 
scales in these cases due to the presence of a Landau pole. 
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In contrast to yLE, the size of of i/le does not significantly affect the maximum achievable 
deviation in Rh^ TT . This is especially true if it is not too large - it contributes exclusively 
through higher dimension operators. These higher dimension operators can cause the con- 
tours of both Ax 2 and Rh^rr to shift, but they move very little with respect to one another. 
So, for a different ule, a somewhat different underlying choice of i/e and yi may be needed 
to achieve a similar effect in Rh^, TT - 

The region of parameter space with both ijle and yLE large can also potentially allow 
for significant enhancement of h — > 77. We return to and elaborate on this point in 
section HVl 



C. Modification of the b Yukawa Coupling 

We now turn to discuss how the measurements of b quarks at the Z pole constrain the quark 
model. To leading order in the couplings we have 

=► SR b = R s b M (l - R s b M )(-1.95gV - 2.7 6g£), 

5A b = O.ffi 6gf- 0.95 5g?, (22) 
=► 5A% = \A s e M 5A b . 

The factor of (1 — R b M ) in the top equation comes from the modification of T(Z — > hadrons) 
via the change in YiZ — > bb). The experimental results and SM predictions are [15J: 

0.21629 ±0.00066 Rf M = 0.21576 ± 0.00004, 
0.923 ±0.020 Af M = 0.9348 ± 0.0001, (23) 

0.0992 ±0.0016 Afi b = 0.1034 ± 0.0007. 

Because Apff deviates from the Standard Model expectation by 2.6<r (leading to Axg M = 6.8 
with respect to the best fit (5g b , 5g b ) point), a requirement of a very small A% 2 with respect 
to the global minimum in the (Sg^, 8g b ) plane is very difficult to satisfy. For our purposes, 
we view it as an unreasonable requirement. After all, this model is not designed to rectify 
this apparent discrepancy with the Standard Model (see e.g. [16J). Instead we use the 
following prescription. Neglecting A^'^ the Standard Model fit greatly improves, yielding 
A%g M = 0.95. We therefore require points exhibit small Ax 2 relative to the global minimum 
from R b and A b only - we neglect ^4^^- We have confirmed that points with a big shift in the 
b Yukawa do not produce a significantly worse fit to than the Standard Model. One can 
think of this Ax 2 for these two measurements as approximately representing the goodness of 
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FIG. 5. Contours (solid) of A% 2 = 2.30, 4.61, and 5.99 relative to the best-fit value in the (6g%,6g£) 
plane neglecting A^'^ as a function of yg and yo for M = Mq = Mo = 600 GeV and £/Q£> = 1, 
Vqd = 0. These values of A% 2 correspond to 68.27%, 90% and 95% regions for 2 parameters. 
For reference, the Standard Model {5g^ = 0, 5g^ = or equivalently tjq = 0, Vd = 0) gives 
Axg M = 0.95. Also shown are contours (dashed) of Rh^bb = T(/i — > bb)/FsM(h — > bb). As 
discussed in the text, larger values of i/qd would allow for larger deviations in Rh^bb- 
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fit relative to that of the Standard Model. As we will see, quite large modifications in the b 
Yukawa can be achieved, so even a somewhat more stringent requirement on the A% 2 could 
yield appreciable effects. 

The analog of Fig. [4] for the quark model is shown in Fig. |5j with Mq = Mo = 600 Ge V and 



Vqd = 1, Vqd = (to reduce tension with AT constraints - see section HID). Subject to the 
requirement A% 2 < 5.99, one can achieve Rh^bb ~ 0.55; recalling Fig. [lj such a modification 
would result in an increase of all other branching ratios by a factor of 1.34. A more extreme 
choice, e.g. ygc = |, i/qd = 0, A% 2 < 5.99, allows Rh^bb ~ 0.4. This point would increase 
branching ratios to other final states (including, interestingly, 77) by 50%. As promised, for 
these points there is not a significant degradation in the fit for A^'^: whereas the Standard 
Model expectation differs from the measured value of by 2.6cr, extreme points in Fig. [5] 
(those with A% 2 « 5.99 and Rh^bb ~ 0.55) exhibit a discrepancy at the level of 3.0a. 



D. Constraints from Oblique Corrections 



From comparison of Eq. (11) and Eqs. (18 19), one can see that it is possible to avoid 
constraints on Sg v and 5g A while generating a larger by compensating for small jje/Me 
and Vl/M l with a larger value for y LE . However, there are constraints on the model 
parameters in addition to the non-oblique corrections already discussed - notably the new 
physics can induce corrections to the S and T parameters [T71 Q2] . 

The T parameter corresponds to the effective operator 

O t = |/i f iVf , ( 24 ) 
such that a AT = v 2 ct- This operator receives contributions of parametric size 

CT ~ T» (25 > 

for yi a Yukawa coupling to the exotic states. In particular, there is a contribution that 
goes as Vleqd- To reach the extreme regions of parameter space in which the deviation of 
R h ^ Tr h ^ bb is the greatest while avoiding tension with AT, one can retreat to higher mass 
scales for the vectorlike particles along with a corresponding increase in Dl,q and de,d- 

While these models also generate contributions to the S parameter via the operator 

O s = (hW a h)W« v B^, (26) 

we find that the constraint from S is typically not significant. In fact, due to the correla- 
tion (88%) between these parameters in the precision electroweak fit [TS] . a small positive 
contribution to S typically allows somewhat larger T values. 
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M 

FIG. 6. Curves of AT for models with either a vector like lepton (N c = 1) or a vector like quark 
(N c = 3) and particular choices of the Yukawa couplings between the heavy fermions. In all cases 
we have taken the vectorlike doublet mass equal to the vectorlike singlet mass. 

The experimental bounds on AT, AS at 2a (making a very slight adjustment for rrih = 125.5 
GeV) are [15] 

- 0.08 < AT < 0.23, (27) 
-0.14 < AS < 0.22. (28) 

Values for AT are shown in Fig. [6] as a function of M = M^^q = Me,d for several choices 
of y = Vle,qd and y = yLE,QD, neglecting mixing with SM particles (which is constrained 
to be small). Full expressions for the contribution to T from vectorlike particles can be 
found in [19] and are reproduced for this specific model in Appendix [Bj Formulae for S can 
also be found by suitable modification of the formulae in [19], but the expressions are more 
lengthy. 

For y = y = 1 (iV c = 1), the AT constraint requires M > 300 GeV at 2a (AT = 0.22, 
AS = 0.09) and M> 360 GeV at la (AT = 0.15, AS = 0.06). As can be seen from the 
figure (y = 0, y = 1, N c = 1), one can abate the tension with AT and permit lower mass 
scales by taking y small - this suppresses contributions to ct involving y, and leads to smaller 
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values for ATj^] For instance, one could take y = and increase y to y = § - in the lepton 
case (JV C = 1), this would allow Rh^rr ~ 0.6 consistent with A% 2 < 5.99 (compared with 
Rh^rr ~ 0.7 for yiE = Vle = 1) while increasing the required mass scale to 380 GeV 
(corresponding to a lightest exotic charged lepton mass of 270 GeV). For yiE = 0, i/le = 2, 
Rh^rr ~ 0.5 is allowed, with required masses 700 GeV. The lightest lepton in this case 
will be at 550 GeV, likely out of reach for the LHC. In both of these cases the contribution 
to S is modest, and the AT nearly saturates the relevant bound of AT^ 0.16. 

In the case of the quark model, tension with constraints on AT is increased due to the 



additional color factor N c = 3. Thus, as in section HI C we consider the region of parameter 
space with {/qd ~ 0(1) and yQD small to avoid AT constraints. Furthermore, as we explain 
in section IV, the region yQD ~ i/qd ~ Oil) is less interesting for h — > 77 for the quark 
model than for the lepton model, in part due to the reduced charge. Consequently, increasing 
yQD would mostly serve to increase bounds on Mq, Me>. 

Again, one can take larger values for i/qd to achieve smaller Rh^bb at the price of increasing 
bounds on M. As mentioned earlier, yQD = |, ygc = admits Rh->bb ~ 0.4. The price 
is that the AT constraint requires M <: 680 GeV. Nevertheless, the lightest quark mass 
in this case is 560 GeV, likely discoverable soon. For |/qd = 2, yQD = 0, Rh-^bb can be 
as low as ~ 0.25, but the mass bound increases to Mq = Mr, = M ^ 1250 GeV. Such 
modifications would correspond to increases of other all branching ratios by factors of 1.5 
and 1.75, respectively. 

To reiterate, if one is willing to permit Landau poles at a relatively low scale, one can go 
to large values of i/le,qd and correspondingly larger values of Ml^q, Me,d and yL,Q, Ve,d to 
achieve smaller values for Rh^ T T,h->bb (or, equivalently, larger c^). However, this requires 
the new states to be more massive, making direct discovery more difficult, particularly in 
the case of the vectorlike leptons. 



IV. MODIFICATION OF HIGGS COUPLINGS TO MASSLESS GAUGE BOSONS 
A. Coupling to Photons 

An extension similar to the one considered here was discussed in pjj with an eye toward 
increasing the 77 branching ratio for the Higgs (and more recently in [TOTTfS] ). There, 

3 In fact, for some negative choices of y it is possible to achieve some cancellation between contributions to 

T - we do not concentrate on this region of parameter space as it would lead to a suppression in the 77 

rate, which at present is disfavored by the data. 
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vectorlike leptons were added without substantial (in some cases any) mixing with Standard 
Model leptons. 

In our case, we have shown (by dialing the r-exotic mixing) that sizable deviations to Higgs- 
r effective Yukawa can be achieved for new leptons of essentially any exotic mass. It is 
therefore straightforward to simultaneously enhance the effective Higgs coupling to photons 
by focusing on the low mass region. In effective theory language, we are generating the 
operator 

0, n = tfhF^F^, (29) 
with coefficient c 77 . Parametrically, we expect contributions 

(30 > 

with yi the largest Yukawa coupling in the problem. In particular, since yi and tje are small 
to satisfy constraints on the Z — r — r coupling, we expect dominant contributions from 
Vi = Ule or y LE . In fact, a simple expression for the coefficient of the operator that couples 
a single Higgs boson to photons 

0k ~'~< = le^""^' (31) 

can be derived via the general formula, [9] 

Ch ^ = bl/2 dh^ hg { detM \ M f) > ( 32 ) 

with b\/2 = (4/3)iV c Qj for a Dirac fermion. Using our mass matrix, neglecting the small 
mixing with the Standard Model leptons, we find to leading order in v 2 /M 2 

16 AT rfVLEVLEV 2 /qq\ 

^ = -T n ^'-mmT (33) 

with Ql — 1, N c — 1 for leptons. Maximizing constructive interference with the WW loop 
(destructive interference with the top fermion loop) requires yLE and yiE of same sign and 
large. This further motivates our choice of yLE = Vle = 1 in the plots of Fig. |4j 

Armed with the above effective field theory understanding, we have identified the region 
of parameter space that gives the maximal change in h — > 77. We now proceed to a full 
numerical evaluation of the effects. The modification of the Higgs boson width to photons 
can be written as: 



r(fe -+ 77) 
r S M(^ -» 77) 



f^i(^) + 2 %? 3 (I) 2 A 1/2 (r top ) + E A A 1/2 (ri) * 



9 -WMrw) + 2 -^S(lYA 1/2 (r top ) 



w 



2 



(34) 
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M L = M E [GeV] 

FIG. 7. Contours of Rh^ n = rsMpSry) (dashed) and Ry l ^ rTT (solid) as a function of the mixing 
Ue = —yi and Ml = Me- We fix y^E = yiE = 1, which is allowed by precision electroweak 
measurements for values of M > 300 GeV at 2a. The shaded region is excluded by too large shifts 
5gV , 5<^5 in the coupling of the r to the Z (A^ 2 > 5.99). 



0.109 



E 



2a. 



htjt 



-^1/2(4; 



(35) 



with Tj = 4m 2 /m 2 t , and loop integrals Aj(t) denned in the appendix. In the second line, 
for the particles with all of their mass from EWSB, we have substituted 9h 1 ^ v = '^ bM - — — 
(v = 174 GeV), as well as the values for the SM loops, A w = —8.34, A top = 1.38. 



T(h -> 77)/T SM (^ -> 77) are shown in Fig. [7| for y LE = y LE = 1 as a 



function of M L = M E and y^ = —yE- Also shown are regions excluded by constraints on 
8gX' A - The choice y^ = —yE is motivated by Fig. [4j which indicates that this is the region 
of parameter space that exhibits the largest deviation in Rh_^ TT . As can be seen from the 
figure, it is possible to simultaneously achieve Rh^ T r ~ 0.7 and Ru-^y-y ~ 1.5 consistent with 
experimental constraints. Doing so requires light new states, making it conceivable that they 



16 



may be observed at the LHC - we briefly discuss potential signatures of these light leptonic 
states in section IVl 

Throughout, for simplicity, we have taken the vectorlike doublet and singlet to have a 
common mass scale, Ml = Me = M. It is reasonable to wonder how sensitive our conclusions 
are to this choice. Rh^rr and 5g^' A are (to leading order) functions of yi/Mi (i = L,E), 
so the values they attain are largely unaffected by perturbations from this point - shifts 
in Ml, Me can be compensated by corresponding shifts in Ul^Ue- The situation is slightly 
different for the T parameter: as AT is a measure of mass splitting within the doublet, it 
is more sensitive to Ml than Mg. Thus, it is possible to achieve the same values of AT by 
decreasing Me and increasing Ml by a smaller amount (such that the splitting in the doublet 
decreases). Such movements in parameter space could be used to slightly increase -R/^ 77 
without diminishing the electroweak fit. However, doing so only increases -Rh-> 77 by (9(5%), 
so we consider points with Ml = Me to be appropriately representative of the variations in 
Rh^-yy achievable. 

For similar parameters, the contributions of down-type quarks to Rh^^ n are less important 
by a factor of three. Furthermore, their Mq^ are constrained by collider (and T parameter) 
considerations to be larger than the lepton case. For instance, for Mq = Mo = M and 
I/q-D ~ Uqd ~ C*(l); * ne lightest new state has mass mi ~ M — v. Current bounds 
on vectorlike quarks constrain M J^, 600 GeV in this case [201422] - for such values of M, 
deviations of -Rh_5. 77 from unity are negligible. Consequently, we find that loop contributions 
to h — > 77 from the vectorlike quarks are generally small. A non-trivial enhancement in the 
Higgs branching ratio to photons can still be achieved as a result of the suppression of the 
effective b Yukawa. This approach requires that the Tevatron excess in bb was not due to 
the Higgs boson. Hopefully, searches for hZ with a boosted h — > bb at the LHC will soon 
help shed light on this point. 



B. Coupling to Gluons 



While the vectorlike quarks do not significantly affect h — > 77, they may affect the Higgs 
coupling to gluons through generation of the effective operator 

O gg = tfhG a ^G; u . (36) 

As in the case of photons, we can consider the coefficient of the operator that couples a single 
Higgs boson to gluons 

ol h 

O h99 = ^--^-G^G^, (37) 
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which is generated with coefficient [23] 



(3? 



where by 2 = (2/3) for a Dirac fermion. To leading order in v 2 /M 2 (and neglecting small 
mixing with the Standard Model), the vectorlike quarks generate a coefficient of size 



Ch 



!/!/ 



8 ygpygpy 2 
3 M D M Q 



(39) 



Thus, depending on the relative signs of i/qd and y~QD, the contribution from the new 
vectorlike quarks can interfere either constructively or destructively with the top loop (which 
dominates the Standard Model contribution). In full, the modification of the Higgs boson 
gluon fusion production cross section is 

2 

\ . 2 9hB,s, a / 4 ^ , 

o(gg ->■ h) 



2 9hti 
mt 



p$>u{gg ->■ h) 



2g h 



-A 1/2 {T top ) 



0.512^ 

Bi 



2ghBiBi V 



^1/2(4; 



(40) 



(41) 



Recent studies [HI [2H I2S] have suggested that fits to the data may favor a slight decrease in 
Higgs production via gluon fusion. However, these studies were performed without the recent 
data from ATLAS on WW |26j, which favor a slightly increased rate. In this model, moderate 
enhancement or suppression of gg — > h are both possible. For yg D ~ 0, effects on gg — > h will 
be small, so a large deviation in Rh^bb from unity can be achieved without simultaneously 
affecting the Higgs boson production cross section. Destructive interference with the top 
loop occurs for same sign yq D and ])qd- For yg D ~ yq D ~ 0(1) and Mq = M D > 600 GeV 
(to avoid direct search bounds for the lightest state), one can achieve R gg ~>h = ^(gg —> 
h)/asM{gg — > h) « 0.7 (simultaneous with Rh^bb ~ 0.8). Alternatively, one can generate 
Rgg^h > 1 for opposite sign y QD and y QD . For y QD ~ -y QD ~ 0(1) and M Q = M 



D ~ 



560 GeV (to satisfy bounds on a charge-2/3 quark decaying exclusively to bW [21]), it is 
possible to achieve R gg ->h ~ 1-35. 



V. PHENOMENOLOGY 



We have implemented the lepton model in MadGraph 5 [27] using FeynRules [28]. Production 
rates for the new heavy leptons are shown in Fig. [81 using the CTEQ 611 parton distribution 
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function. For the points we have chosen (ule = Vle = 1,Me = Ml) the masses obey the 
simple relation mi x + 174 GeV = hin = m^ 2 — 174 GeV. Bounds on this model exist from 
LEP near the kinematic limit (m^ > 105 GeV) [29] . 

In general, the lightest charged state, £\, will decay to either rh, tZ or Wv. One potentially 
relevant search is the ATLAS slepton/chargino dilepton search [30]. But while it overlaps the 
final state, it is not yet sensitive, although future searches might be at low masses. Notable 
is the rh final state |31j . which is particularly important in this region of parameter space. 
The presence of both a vectorlike L, L and E, E allows deviation from the characteristic 
1:1:2 ratio of final states found for models with only a singlet vectorlike partner (i.e. a 
"Littlest Higgs"-like model). For example, for points with i/le = Vle = 1 (as shown in the 
plot) and He = —Vl large (so as to give a large suppression in Rh^ TT ), the branching ratio 
BR(£\ — > hr) is near 50%. Of course, for this final state to be relevant, the lepton must 
be heavy enough to evade the phase space suppression, which limits its production cross 
section. Nevertheless, we find the BR is nearly 50% already at mt x = 175 GeV. With high 
luminosity, one might even explore the possibility of utilizing h — > 77 decays along the lines 
of [32], particularly given the possible enhancement of the 77 rate. The N state will almost 
exclusively decay to £iW, raising the possibility of multi-lepton cascades. 

In the case of vectorlike quarks, the phenomenology is similar to well-explored heavy quark 
models - production is dominated by QCD processes, and cross sections can be determined 
as a function of mass [331134]. We review some of the relevant limits here (see also [35]). CMS 
has searched for a B' via B' — > bZ. For a B' with a 100% BR to bZ, the limit is > 550 
GeV p3]. A similar search from ATLAS, but using less data, sets a limit > 400 GeV 
[36] . Using 4.9 fb _1 of data the CMS collaboration excludes a B' decaying to tW with 100% 
BR below 611 GeV at 95% confidence [20]. 

Searches for T — > tZ by CMS |37J exclude a top partner with 100% BR to tZ at 95% CL. 
Searches for bW final states (with two leptonically decaying W's) |21j exclude masses up 
to 557 GeV. A search in the semileptonic final state by ATLAS [3H] has a more limited 
reach of 480 GeV. A combination and reinterpretation of these searches can be used to 
bound a T with non-trivial branching ratios to all of tZ, bW and th |22j. Some degradation 
of the above limits exists, perhaps by up to 100 GeV or so. A dedicated search for th 
would improve the situation. In our case, depending on the exact implementation of the top 
sector, the T will likely first cascade to a B\W, followed by a further decay to SM fermions 
and electroweak bosons. This opens the possibility of, e.g., WWZZbb final states. For 
Vqd ~ —Uqd ~ (such that gluon fusion gg — > h is enhanced), T can be the lightest new 
state and will consequently decay to tZ, bW and th - neglecting mixing in the top sector, T 
decays exclusively to bW and the bound of ttit > 557 GeV applies. More detailed collider 
studies are left for further work. 
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FIG. 8. Production cross sections for the heavy lepton states for y^E = Vle = 1 and Mi = Me 
at the LHC with yfs = 8 TeV (top) and 14 TeV (bottom) as a function of the lightest exotic 
lepton mass. Along the top axis the enhancement in the 77 branching ratio with respect to 
the Standard Model value is shown. Note for yiE = yiE = 1; constraints on (AT, AS) require 
Ml = Me > 300 GeV, or > 126 GeV. These plots are valid for any values of vliVe consistent 
with constraints on 5g^ , ^9r- 
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FIG. 9. The dashed lines indicate the one-loop running of the gauge couplings in the extension 
considered here with three pairs of vectorlike leptons. The solid lines indicate the running in the 
Standard Model. 



VI. COMMENTS ON UV COMPLETIONS 



One possibility is that the model presented here might be embedded in a (perhaps somewhat 
split [391FTT] ) SUSY scenario, where the scalars superpartners are sufficiently heavy that they 
do not affect the phenomenology discussed here. In this case, approximate gauge coupling 
unification is maintained only via the introduction of complete £77(5) multiplets. With a 
single 5 + 5 and 10 + 10 perturbative gauge coupling can be maintained. So, introduction of 
the vectorlike leptons would also motivate the presence of the vectorlike quarks - this would 
be good news for LHC phenomenology. In our analyses in the previous sections, we have 
not adhered too strictly to this motivation, as it would imply relationships between Yukawa 
couplings in the (exotic) lepton and quark sectors that we have not imposed. The C(TeV) 
masses considered here could be explained by whatever mechanism is responsible for the 
solution to the \i term (e.g. Giudice-Masiero or NMSSM like physics). 

It is amusing to note the possibility of a less conventional unification story. Suppose only 
vectorlike leptons are added (without the corresponding quarks). With three pairs (i.e., 
3 x (L,L,E,E)), as might be expected if there is a vectorlike partner for each generation, 
there is adequate unification, see Fig. [9j Unification can be assessed by examining the ratio 
R = (b 3 — b 2 )/(b2 — bi), where 6j denotes the /3-function for the gauge group i. Under the 
assumption of unification, at one loop R = (a^ 1 — &2 )/{ a 2 ~~ a \ 1 )- Experimentally, the 
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RHS evaluated at Mz yields 0.718 ± 0.003. In the Standard Model, R = .528, indicating an 
exceedingly poor fit for unification. Augmenting by three vectorlike pairs, we find R = .707, 
which while could plausibly be accounted for by additional corrections at the GUT scale. 
The scale of unification is quite low Mjj = 2.4 x 10 10 GeV, so additional GUT model building 
(likely with some type of extra-dimensional unification) would be necessary to avoid too large 
dimension six proton decay. However, it should be noted that this approach to unification 
creates new doublet-triplet splitting problems. In addition, for larger Yukawa couplings 
(where the largest effects on the Higgs branching ratios are achieved), there is a danger of 
introducing vacuum instability with decay times shorter than the observed lifetime of the 
universe [TT] . 



VII. CONCLUSIONS 



We have explored the possibility that vectorlike fermions can mix with the Standard Model, 
giving substantial modification to the Higgs boson properties. To summarize the types of 
deviations possible, we have constructed Table|TJ which illustrates the extent to which R h ^ TT , 
Rh^bb and Rh-^y-y can be modified. While the points in this table represent some of the more 
extreme cases (even allowing a total enhancement of .R/i_> 77 ~ 2), it should be noted that 
there are any number of effects in play here that can work in concert, all quite plausibly 
present if there is new vector like matter with significant Yukawa couplings at the TeV 
scale. 

For example, the observed a x BR for photons can be affected in multiple ways: vectorlike 
quarks might enhance the gluon production cross section, the presence of vectorlike leptons 
might increase T(h — > 77), and mixing with the b quark (r lepton) could reduce T(h — > 
bb) (r(/i — > rf)). These effects can be sizeable: we have found realizations where Th-^rf 
is suppressed by a factor of 2, where T h ^ bb - is suppressed by a factor of 4, where r^ 77 
is enhanced by 50% and where cr gg ^h is enhanced by 35%, all consistent with precision 
constraints. Some of these effects push against the same experimental limits, so while some 
variations can be thought of independently, others (such as simultaneous enhancements to 
Ogg^h an d r^^p,) cannot be. Nonetheless, one can find points where the overall inclusive 
signal of h — > 77 is enhanced by a factor of 2. 

It is important to emphasize that all these effects are naturally present in models with 
vectorlike fermions. One need not focus on pushing every mode to the limit to produce an 
interesting effect. Even if individually modest, multiple contributions can combine to give 
an interesting effect. A combination of ~ 10 — 20% effects to yield a 50% modification is 
easily plausible without pushing the limits of existing constraints. 
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fa) R 



h — }TT 



0.53 



M L = M E [GeV] 


VLE 


VLE 




VE 


AT 


AS 


A X 2 r 


Rfl^fTT 


Rh—^yy 


Masses [GeV] 


700 


2.0 


0.0 


0.105 


-0.105 


0.16 


0.01 


5.98 


0.53 


1.00 


547, 700, 895 


(b) Rh^TT = 0.74 and Rh^yy = 1.50 


M L = M E [GeV] 


VLE 


VLE 


2/L 


VE 


AT 


A 5 


Ax? 


Rfl-^TT 


Rh^f-yy 


Masses [GeV] 


300 


1.0 


1.0 


0.032 


-0.028 


0.22 


0.09 


5.90 


0.74 


1.50 


126, 300, 474 



(c) R, 



h^bb 



0.25 



M Q = M D [GeV] 


VQD 


VQD 


VQ 


VD 


AT 


AS 


AX^ 


Rh-^bb 


non— b 


Masses [GeV] 


1250 


2.0 


0.0 


0.57 


-0.27 


0.16 


0.02 


5.93 


0.25 


1.74 


1091, 1250, 1437 



(d) 



BR(fe-»77) 
BR SM (^77) 



1.99 



M Q = M D [GeV] 


VQD 






VD 


AT 


AS 


AX^ 


Rh^bb 


non- 6 


Masses [GeV] 


1500 


1.5 


0.0 


0.71 


-0.31 


0.04 


0.01 


5.86 


0.39 


1.53 


1378, 1500, 1638 


M L = M E [GeV] 


VLE 






VE 


AT 


AS 


Ax? 


Rh-^rTT 


n h^yy 


Masses [GeV] 


350 


1.0 


1.0 


0.0 


0.0 


0.16 


0.06 


Ax| M 


1.00 


1.30 


176, 350, 524 



TABLE I. Benchmark points exhibiting (a) large deviation in R^ TT from unity (lepton sector), 

(b) moderate deviation in Rh^ TT from unity and maximal value of R^^ from the lepton sector, 

(c) large deviation in Rh^bb from unity (quark sector) and (d) large enhancement of h — > 77 due to 
combination of additional loop contributions from vectorlike leptons and suppression of b Yukawa 
from mixing with vectorlike quarks. r^_ b = BR(h — > X) / 'BR$yi(h — > X) for X / bb as a result 
of the decrease in the h —> bb width. Note that, for points (a) and (c), we have taken M as small 
as possible consistent with EWPT constraints - for larger values of M, the same Rh^ T T,h^bb could 
be achieved with smaller values of AT, AS. For point (d), -R^-m denotes the enhancement of 
h — > 77 from the lepton sector only. 
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All but the most extreme deviations in the fermion couplings discussed here will be challeng- 
ing to definitively confirm at the LHC [42J , but would be easily measured at a linear collider. 
For a recent overview, see |43j . While deviations in the Higgs couplings motivate this model, 
potentially just as interesting are the direct searches for the new states introduced here. The 
largest modifications to the 77 rate are only possible for for light leptons with masses, very 
likely accessible soon at the LHC. 

A similar Lagrangian was considered in [H] in an attempt to explain the PAMELA positron 
excess. There, a single vectorlike lepton was introduced (rather than both E and L). (In that 
case, all corrections are proportional to y®, and will not be numerically significant.) However, 
just as in that case, if the Dark Matter (here left unspecified) couples to the vectorlike leptons 
introduced here, it could lead to relatively leptophilic annihilation channels, with important 
implications for indirect detection. 

The discovery of the Higgs boson is a watershed moment in particle physics. If new physics 
exists at the weak scale, it may leave imprints on the measurable properties of the Higgs 
boson. What is remarkable about a 125 GeV Higgs is that, with many observables, its 
properties are so sensitive to many scenarios of BSM physics. The presence of new vectorlike 
fermions is intriguing, because a complete generation can influence independently r^ 77 , 
gg-fusion, Th^ Tf and T h ^ b i. Should future data hold up an anomaly in these quantities, it 
may point to new matter just around the corner. 

Note Added: While this paper was nearing completion, some works pTJI4T2] ) appeared that 
discussed the enhancement to 77 due to vectorlike fermions. Also, [45] noted the possibility 
of suppressing the r Yukawa via a mixing with vectorlike leptons in a similar model with an 
additional singlet scalar. 
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Appendix A: Higgs Loop Functions 



For completeness, we reproduce the loop functions for the generation of the h — > 77 process 
used in Sec. [TV] 

-VM=2r 2 (i+(i-l)/(i)), (Al) 
^W = -^(^ + |+3(|-l)/(i)), (A2) 

,2 



f(x) = Arcsin 2 A/x. (A3) 

The form of f(x) is valid for rrih < 2mj, where rrii is the mass of the particle in the loop, as 
for the particles considered in this paper. 



Appendix B: Functions for the T parameter 

Neglecting mixing with the Standard Model particles (which is in any case constrained to 
be small), we find a contribution to the isospin breaking T parameter as [H]. In this limit, 
the mixing of the charged vectorlike leptons is described by the 2x2 mass matrix 

- £i,,,ss ( E<- I ) I ME VLEV I I ~ 1 +I1.C. (Bl) 



y LE v M L 

which can be diagonalized to yield mass eigenstates 

E c 








and _ = U* _ , (B2) 



where M.d = diag(rni,W2) = UJ\AV\ The spectrum also contains of a vectorlike neutral 
state with mass Ml. Modifying the expression from [19J, we find the contribution to the T 
parameter 

AT = 1R N 2 2 \ E fd^l 2 + \ U *\ 2 ) d + (Vi'Vi) + 2Rfi(VbC^)0_ ( yi ,y L ) ) 
wws w c w ^ . \ / 

\V 12 V 22 \ 2 + \U 12 U 22 \ 2 ) 9 + ( yi ,y 2 ) + 2Re(V 12 V* 2 U* 12 U 22 ))9_ (y u y 2 ) ) J (B3) 

where yi = mf/m 2 z , N c is the number of colors and 

9+ {vu V2) =yi+V2 log— (B4) 

yi - 2/2 2/2 
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e-(vi, a.) = ( VA± ^- log - - 2V (B5) 

Completely analogous expressions hold for the quark model. 
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